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Context. Many galactic B[e] stars suffer from improper distance determinations, which make it difficult to distinguish between a pre- 
and post-main sequence evolutionary phase on the basis of luminosity arguments. In addition, these stars have opaque circumstellar 
material, obscuring the central star, so that no detailed surface abundance studies can be performed. 
' Aims. Instead of studying the surface abundances as a tracer of the evolutionary phase, we propose a different indicator for the 

supergiant status of a B[e] star, based on the enrichment of its circumstellar matter by 13 C, and detectable via its 13 CO band emission 
Q in the K band spectra. 

Methods. Based on stellar evolution models, we calculate the variation of the 12 C/ 13 C isotopic surface abundance ratio during the 
evolution of non-rotating stars with different initial masses. For different values of the 12 C/' 3 C ratio we then compute synthetic first- 
overtone vibration-rotational band spectra from both the 12 CO and 13 CO molecule at different spectral resolutions. We further discuss 
the influence of stellar rotation on the variation of the surface 12 C/ 13 C ratio and on the possibility of 13 CO band detection. 
Results. The surface 12 C/ 13 C isotope ratio is found to decrease strongly during the post-main sequence evolution of non-rotating stars, 
from its interstellar value of about 70 to a value of about 15-20 for stars with initial masses higher than 7 M Q , and to a value of less 
than 5 for stars with initial masses higher than 25 M Q . We find that detectable l3 CO band head emission is produced for isotope ratios 
12 C/ 13 C S: 20, and can most easily be detected with a spectral resolution of R ~ 1500 . . . 3000. For the rotating stellar models, the 
drop in 12 C/ I3 C already occurs for all stars with M m i, 9 M Q during the main-sequence evolution. The detection of 13 CO band head 
emission in such mid-resolution K band spectra of a B[e] star thus favours an evolved rather than a young nature of the object. 
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\ 1. Introduction their winds comes from polarimetric observations, e.g., by 

Magalhaes d!992l l. Magalhaes et al. (|2006), and Melgarejo et al. 

O ■ The classification of B[e] stars according to their evolutionary (2001 3 for the Magellanic Cloud stars, and from, e.g., Zickgraf 

L; ; phase is a long-standing problem. Much effort has been under- & Sc hulte-Ladbeck ( 1983 for some of the galactic candidates. 

. . _h . taken since Lamers et al. 0558) published the first detailed clas- Detection of mo i eC ular emission bands of TiO (Zickgraf et 

^ ■ sification criteria sorting the B[e] stars into the classes of su- a i.[j989]l and CO (McGregor et al. [Mil DSUl Morris et al. 

»H i pergiants, Herbig stars, compact planetary nebulae, and sym- m6) in the spectra of several B[e] supergiants connrms the 

<S? , biotics. Nevertheless, the remaining group of unclassified B[e] presence of high-density and cool material in the vicinity of 

stars, which includes about half of all known B [e] stars, has since th^g luminous stars 
then gained even more members. 

While many targets within the group of unclassified B[e] An additional characteristic of B[e] supergiants isjheir 
stars have not been studied yet in detail, there also exists a num- stron g lnfrared ( IR > excess emission (e.g., Zickgraf et al.[T98l, 
ber of stars showing indications of both a young (Herbig star) as indicating that these luminous objects are surrounded by hot 
well as an evolved (supergiant star) nature. And one of the rea- or warm dust ' whose most probable location is within a mas- 
sons for their uncertain classification is their still rather poorly Slve circumstellar disk. Recent estimates, based on observa- 
known distance, and hence luminosity (see TableED- For a tlons wlth the S P ltzer S P ace Telescope Infrared Spectrograph for 
proper assignment of an evolutionary phase to these stars it is, the Lar 8 e Magellanic Cloud B[e] supergiant R 126 by Kastner 
therefore, necessary to find other characteristics that are directly et aL I 2 ™' resulted ln a total dust mass Wlthln the dlsk of 
linked to the evolutionary phase of a star as additional classifica- ~ 3 x 10 M . 

tion criteria. The existence of dusty disks is also one of the major char- 

The best studied B[e] supergiant sample is located in the acteristics of Herbig Ae/Be stars (see, e.g., the recent reviews 

Magellanic Clouds. Their main characteristic is the so-called by Waters 2006 and Waters & Waelkens 1998 and references 

hybrid character of their spectra, showing indications for both therein). And for many young stellar objects the observations of 

a typical line-driven wind in polar direction, and a slow and the CO band heads (Hanson et al. 119971 Bik & Thi 120041 Thi 

low-ionized but high-density wind in the equatorial direction et al. 2005] Bik et al. 120061 1. although not present in all objects, 

(Zickgraf et al. l 19851 >. Further evidence for non-sphericity of often resulted in detailed information about the disk structure 



2 



M. Kraus: The pre- versus post-main sequence evolutionary phase of B[e] stars 



and kinematics of the inner disk. A high-density circumstellar 
disk with molecules and dust is thus a common feature of both 
the Herbig stars and the B[e] supergiants. The only, and conse- 
quently most important difference is the fact that for the Herbig 
stars the disk material is the remnant of the star formation pro- 
cess. 

Since the progenitors of B[e] supergiants must have been 
massive O-type stars, this close-by dust cannot be pre-main se- 
quence in origin. Any remaining pre-main sequence circumstel- 
lar material will have been swept away by the strong radia- 
tion pressure of their stellar winds during their main-sequence 
life-time. The disks around the B[e] supergiants must thus have 
formed from the wind material itselfl due to some enhanced 
mass loss or some violent high-density mass ejection, predom- 
inantly within the equatorial regions. Detailed investigations of 
the UV metallic lines for some Magellanic Cloud B[e] super- 
giants seen edge-on have revealed terminal disk wind velocities 
of 60-80kms~ 1 only (Zickgraf et al. 1996). They are thus at 
least a factor of ten lower than the terminal velocities in the polar 
winds, which show properties similar to those of classical B-type 
star, i.e., line-driven winds (e.g., Zickgraf et al. 1985). In com- 
bination with the enhanced mass fluxes found for the equatorial 
wind, B[e] supergiants show a density contrast of 100-1000 be- 
tween the equatorial (i.e. disk) wind and the normal (i.e. polar) 
wind (Zickgraf et al. 119891 ). Reasonable attempts at explaining 
these equatorial flows are provided by the rotationally induced 
bi-stability mechanism (Lamers & Pauldrach 1991; Pelupessy 
et al. 120001 Cure |2004l Cure et al. 120051 Even though the for- 
mation mechanism is not fully understood and solved yet (see, 
e.g., Ignace & Gayley 2005, Kraus & Miroshnichenko 2006), 
the disk-forming winds of the B[e] supergiants must be dense 
and cool enough to allow for efficient molecule and dust forma- 
tion. 

The most stable molecule is CO, and during their IR spectro- 
scopic studies of early-type emission line stars, McGregor et al. 
( 1988 ) detected CO band head emission from several targets that 
later grouped in the unclassified B[e] stars or as B[e] supergiant 
candidates (see Table lA.lt . 

What remains is the question of where in the disk, i.e., at 
what distance from the B[e] supergiant star, the CO band emis- 
sion originates. Recently, Kraus & Borges Fernandes (2005|) and 
Kraus et al. (2006, 2007 ) have shown that the disks around B[e] 
supergiants must be predominantly neutral in hydrogen close to 
the stellar surface in order to reproduce the observed strong [Oi] 
emission lines. The temperature in the [Oi] line forming regions 
ranges from ~ 8000 K to ~ 6000 K, while the dust must be lo- 
cated further away from the central star, where the disk tempera- 
tures dropped below the dust evaporation temperature of 1500 K. 
CO has a dissociation temperature of about 5000 K. Theoretical 
CO band calculations have shown that the formation of pro- 
nounced band head structures requires a minimum CO tempera- 
ture of about 2000 K (see, e.g., Kraus [1997l and Sect.lXTlbelow). 
The CO bands, therefore, probe the disk area between the [Oi] 
and the dust emission regions. In terms of distance, this temper- 
ature argument places the CO formation region beyond the [Oi] 
saturation, which usually happens at distances of around 100- 
300 R, (see Kraus et al. l2006ll2007l l. 



1 This argument is only strictly valid if we assume that the B[e] phe- 
nomenon is a special phase (for stars with perhaps very specific initial 
conditions) during single star evolution, which we will assume through- 
out this paper. If the B[e] phenomenon for these supergiants is caused 
by some binary interaction or binary merger process, the situation might 
be completely different. 



Since the disks around B[e] supergiants are formed from 
wind material, we can expect the disk material to mirror the sur- 
face composition of the star at the time of ejection. Depending 
on the evolutionary phase of the star at the time of matter ejec- 
tion, the disk material already should be enriched in chemically 
processed material that has been brought to the stellar surface 
via rotation and mixing processes. Searching for such enhanced 
abundances within the disk material of B[e] stars might thus pro- 
vide us with the information necessary to discriminate between 
a pre- and post-main sequence evolutionary phase of the under- 
lying star. 

One of the elements produced during the evolution of mas- 
sive stars is the carbon isotope 13 C. Its ratio with the main iso- 
tope 12 C changes drastically during stellar evolution, making 13 C 
an ideal tracer for post-main sequence evolutionary phases. We 
already know of the existence of 12 CO band emission from the 
disks of B[e] supergiants, but none of the targets has been re- 
ported to show 13 CO emission so far. The reason for this is 
certainly the fact that either the spectral resolution used was 
generally too low (R £ 500, see Sect.O, or, when using high- 
resolution spectrographs, the spectral region covered by the ob- 
servations where restricted to the region around the first 12 CO 
bandhead only, which does not include the range of 13 CO emis- 
sion. This paper, therefore, investigates the possible appearance 
of 13 CO band head emission during the evolution of OB stars, 
based on theoretical model calculations of the first-overtone 
bands of I2 CO and 13 CO in the K band spectra. 

The paper is structured as follows: In Sect.|2] we study the 
variation of the 12 C/ 13 C surface abundance ratio during the evo- 
lution of predominantly massive stars. Next, we calculate sim- 
ple synthetic CO first-overtone band spectra in Sect.[3]of the CO 
molecul^l and discuss the appearance of the band heads of the 
I3 CO molecule as a function of 12 C/ 13 C ratio as well as of spec- 
tral resolution. In Sect.|4]we discuss the reliability of our calcu- 
lations and applicability to the B[e] supergiants, before we sum- 
marise our results in Sect. [5] 



2. Variation of the C surface abundance during 
stellar evolution of massive stars 

During the evolution of massive stars, heavy elements produced 
during the different burning processes are brought to the stellar 
surface via rotation and mixing processes (see, e.g., the reviews 
by Pinsonneault 119971 and Maeder & Meynet 2000). The sur- 
face abundances of individual elements will thus change with 
time. Due to the (high) mass loss of massive stars via line-driven 
winds, these changes in surface abundances will thus translate 
into abundance changes in the wind, which then clearly shows 
the presence of chemically enriched, processed material. Studies 
of the chemical composition of the wind (or more generally: of 
the circumstellar) material of post-main sequence evolutionary 
phases of massive stars are, therefore, ideal to determine the evo- 
lutionary phase of the central star. 

To study the changes in surface abundances (and hence of 
the wind abundances), we use the grid of Schaller et al. (119921 ) 
for non-rotating stars at solar metallicity. In Fig. Q] we plot the 
changes in surface abundances along the evolutionary tracks (in- 
dicated by the grid numbers) of 12 C, 13 C, 14 N, and le O with re- 
spect to the surface hydrogen abundance at the same grid point. 
The ratios have been calculated from the mass fractions provided 



2 Note that by writing CO without indication of the isotope, we al- 
ways refer to the main isotope, i.e., the I2 CO molecule. 
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Fig. 1. Changes in 12 C, 13 C, 14 N, and 16 surface abundances with respect to the hydrogen surface abundance along the evolutionary 
tracks of non-rotating stars, as indicated by the grid points. The ratios were calculated from the data of Schaller et al. (j 1 992) for 
solar metallicity. The approximate end of the core H burning and the onset of the core He burning phases are indicated by the two 
vertical dotted lines in each panel. 
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Fig. 2. Change of the 12 C/ 13 C surface abundance ratio during stellar evolution of non-rotating massive stars at solar metallicity. The 
ratios are calculated from the tables of Schaller et al. Q992 ). The dotted lines have the same meaning as in Fig.[T] 



in the data files of Schaller et al. (1992). The grid points that in- 
dicate the approximate end of the core H burning and the onset 
of the core He burning phases are indicated by the two vertical 
dotted lines in each panel. The abundance curves for stars with 



initial masses M; n i > 40 M Q stop around grid point ~ 30 (i.e. be- 
fore the end of core He burning) since these stars have lost their 
complete outer hydrogen shell via strong mass loss. 
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The change in chemical surface composition during the 
main-sequence lifetime is quite obvious for the most massive 
stars, while the surface abundances of the less massive stars re- 
main at their initial values until the onset of core He burning or 
even longer. The elements with clearly increasing abundances 
during the evolution are 14 N and I3 C, while at the same time the 
abundances of ls O and 12 C slightly decrease. The nitrogen en- 
richment of the surface (and hence the wind) material is usually 
used as an indicator to trace evolved evolutionary phases of mas- 
sive stars, as for the classical OB supergiants (see, e.g., Crowther 
et al Ull Hillier et al. [20031 Evans et al. 120041 Searleetal. 
2008) or the Luminous Blue Variables (e.g., Lamers et al. 2001). 

The direct study and comparison of the surface abundances 
of the B[e] stars would be an ideal tool for the determination 
of their evolutionary phases. Unfortunately, due to their usually 
high density circumstellar material (caused by high mass loss 
in the case of B[e] supergiants), their spectra are often crowded 
with pure emission lines, and the photosphere of most targets is 
not visible at all. A direct way to measure their surface abun- 
dances is thus not possible. The only way to draw conclusions 
about their surface composition is to study the chemical abun- 
dances from their emission lines formed within their close-by 
circumstellar material. However, to determine the CNO abun- 
dances of B[e] stars from their emission lines, it is necessary to 
have a detailed description of the exact geometry, the inclination 
angle of the system, and the latitude dependent ionization struc- 
ture. This information is usually not available and can hardly be 
derived from the modeling of the emission lines alone. We there- 
fore need another, more reliable indicator for the evolutionary 
state of an individual object. 

One important indicator of post-main sequence evolutionary 
phases is the surface abundance of the isotope 13 C, or more pre- 
cisely, the change in the isotopic ratio 12 C/ 13 C. The strong in- 
crease in 13 C in combination with the (slight) decrease in 12 C 
surface abundances (see Fig.[TJ results in a rather steep decrease 
of the 12 C/ 13 C ratio from its initial interstellar value of ~ 70 to 
values < 20 for stars with initial masses ^ 40 M Q , and for stars 
with initial masses ^ 60 M even to 12 C/ 13 C 5= 5. The variation 
in the 12 C/ 13 C surface abundance ratio during the evolution of 
massive stars at different metallicity is depicted in Fig. [2] 

With an abundance ratio lower than about 20, the 13 C iso- 
tope should be easily detectable. However, our intension is not to 
search for individual emission lines of the 13 C atom, but rather to 
concentrate on the molecular emission of CO and its 13 CO iso- 
tope. Besides, our investigations will be restricted to the 12 CO 
and 13 CO, only, since the abundances in 17 and ls O compared 
to the main isotope 16 are generally extremely low, even during 
the late evolutionary phases of massive stars. Contributions from 
the C n O and C 18 isotopes to the total CO band structures are 
thus negligible. 

3. CO band formation 

CO first-overtone band emission is a well known feature arising 
especially from the rotating accretion disks around young stellar 
objects (YSO), and the physics of CO band formation in such 
rotating pre-main sequence disks has already been discussed 
widely in the literature (see, e.g., Scoville et al. 1980; Carr 1989; 
Najita et al. [19961 Kraus |2U00l Kraus et al. 120001 . The main fo- 
cus in YSO research is therefore especially the structure of the 
2 — > band head, which contains the complete velocity infor- 
mation, and from which the rotation velocity, disk inclination, 
as well as some information concerning the density and temper- 
ature within the CO forming disk region can be derived (e.g., 



Carr et al. [19931 C handler et al. IT9931 Carr [19951 N ajita et al. 
[19961 Krau s |20001 Kraus et al. 120001 Bik & Thi |2004l Berthoud 
et al. l2007l 

Such high-resolution 2 — > band head observations also 
would be very important for B [e] stars in studying the kinematics 
of the CO gas, especially for the discrimination between whether 
the disks around B[e] supergiants are Keplerian rotating (i.e., 
stable, long-lived disks), or outflowing disk-forming winds (see, 
e.g., Porter 2003; Kraus et al. 2007). For our purpose related to 
the B[e] supergiants' evolutionary phase, however, we need to 
study a much larger (in wavelength) portion of the complete CO 
first-overtone band structures, because we want to find out (i) 
where in the K band spectra of B[e] supergiants we can expect 
the 13 CO band heads to show up, (ii) during which evolutionary 
phases observable 13 CO band head emission might be generated, 
and (iii) what the most suitable spectral resolution is to easily 
detect these 13 CO band heads. 

3.1. Synthetic band spectra of 12 CO 

For the calculation of the CO band spectra, we follow the de- 
scription of Kraus et al. (2000). Due to the lack of knowledge 
of the proper shape and structure of the B[e] supergiant stars' 
disks, we keep the calculations simple, assuming the CO gas to 
be isothermal. We further assume that the CO gas is in local 
thermodynamic equilibrium. Then, the emission at frequency v 
is simply given by 

I v = B v (T C0 )(l-e- T *) (1) 
with the optical depth 

Tv = K y n C oNn co . (2) 

Here, k v ii co is the absorption coefficient per 12 CO molecule, and 
Nn co is the column density of the 12 CO gas for which we adopt a 
value of 10 17 cirT 2 in order to keep the emission optically thin at 
all wavelengths over the complete band spectrum. We take into 
account a total of 10 vibration bands, and within each vibration 
band a total of 140 rotation levels. The so-called Dunham coeffi- 
cients to calculate the energy levels are taken from Farrenq et al. 
dl991| ), and the Einstein coefficients for spontaneous emission 
from Chandra et al. (II 9961 ). 

To calculate k v n co , we need to specify the profile function 
of the individual vibration-rotational lines for which we use, for 
simplicity, a Gaussian profile. The width of the Gaussian we de- 
scribe by a fixed (internal) CO velocity of 30 km s _1 . The choice 
of this velocity is not crucial, unless the internal velocity exceeds 
the one given by the spectral resolution. For some first test cal- 
culations, we smear the CO band spectra to an arbitrary spectral 
resolution of R = 1500, corresponding to a velocity resolution 
of about 200 km s -1 . 

To study the variation in the individual band head peak in- 
tensities, we calculate the structure of the first-overtone bands 
for different CO temperatures, ranging from 1000K to 5000 K. 
The results are shown in Fig. [3] with the positions of the band 
edge^l indicated by the ticks in the top panel. The formation of 
clear band heads starts for temperatures 7co ~ 2000 K, while for 
CO temperatures below 2000 K no individual and clearly iden- 
tifiable band heads emerge. For Tco — 1000 K only the levels 
within the first vibrational band are excited, giving rise to the 

3 For clarification, we use the terminus band head when talking about 
the intensity peaks, while the transition of lowest wavelength in each 
band we call the band edge. 
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Fig. 3. Influence of the temperature on the synthetic 12 CO band 
emission spectra. Clear band head structures form for Tco > 
2000 K. 



broad, double emission feature. The higher the temperature, the 
more vibrational bands can be excited and the more individual 
band heads appear. A further effect of incresing CO temperature 
is the strong increase in band head intensities, with the 3 — > 1 
as the clearly dominating one, while the 2 — > band head be- 
comes depressed with increasing temperature due to depopula- 
tion into the higher vibrational levels. Due to this strong intensity 
increase with CO temperature, most of the observable band head 
emission originates from the hottest CO gas, i.e., from the disk 
regions closest to the star. 



3.2. The appearance of the 13 CO bands 

The slightly higher mass of the 13 C isotope has two conse- 
quences for the energy levels within the diatomic molecule: (i) 
A higher mass results in a higher moment of inertia, reducing 
the energy of the individual rotational levels within each vibra- 
tional band, (ii) Since the frequencies of the vibrations depend 
on the reduced mass only, the vibrational energy levels within 
the more massive I3 CO isotope will all be reduced compared to 
those within the 12 CO molecule. Consequently, the net effect of 
the higher mass of the isotope is a shift of the coupled vibration- 
rotational transitions, and hence of the complete band spectrum, 
to longer wavelengths. Since the wavelength shift (i.e., the re- 
duction in energy) is especially strong for the vibrational transi- 
tions, a clear wavelength separation appears between the bands 
of 12 CO and 13 CO. The lowest band edge of the first-overtone 
bands arises at 2. 3448 //m (2 — > 0), i.e., at a wavelength be- 
tween the second and third band edge of the 12 CO bands. This 
fact is of major importance for the identification of the I3 CO 
band heads in the total spectrum. 

Since the main effect of the different masses is a shift of the 
individual vibration-rotational lines to higher wavelengths, the 
bands of the 13 CO molecule will, in general, have a very similar 
shape to those of I2 CO. To compute the total band spectrum con- 
sisting of the 12 CO and I3 CO transitions, we again calculate the 
intensity, I v , at each wavelength with Eq. ([TJ, but for the optical 
depth, t v , we now need to include the contribution of the 13 CO 
molecule into Eq. ©, giving 



2 CO jVl2 CO • 



3 CO jVl3 CO . 



(3) 



where M3 C0 and k v n co are the column density and the absorp- 
tion coefficient per 13 CO molecule, respectively. 

It is convenient to parametrize the column density of the iso- 
tope, -/Vn co , m terms of Nn co and the carbon isotope abundance 
ratio, ni2 C /ni3 C , for which we assume that it equally translates 
into the CO isotope ratio, so that we can write 



Nuco = Nn co = Nuqo = 



Nl2 CQ 



nn C o 



Nn CO 

nn c /nn c 



(4) 



with the carbon isotope abundance ratio, nnc/nn c , as the only 
free parameter. From this relation it is already obvious that the 
emergence of the 13 CO bands strongly depends on the carbon 
isotopic ratio, i.e., the lower this ratio, the higher the probability 
of 13 CO band appearance. 

In our test calculations we use the following values for the 
«i2 C /«i3 C abundance ratio: 70, 15, 10, and 5. While the first 
value approximately represents the interstellar (i.e., initial) ra- 
tio, the latter appear during the more interesting, later evolution- 
ary stages of stars with different initial masses (see Fig.[2j». The 
emerging total spectra for these four ratios and for different spec- 
tral resolutions are displayed in Fig.|H with the wavelengths of 
both, the I2 CO and 13 CO band edges, indicated in the top panels. 

For an interstellar carbon isotopic ratio of 70, the 13 CO bands 
remain undetectable. The first indications for 13 CO band heads 
peaking out of the total CO spectrum are for an isotopic ratio of 
15, while for even lower ratios, the individual 13 CO band heads 
are clearly present. However, the appearance of these band heads 
strongly depends on the chosen spectral resolution. To show the 
effect of spectral resolution, we calculated synthetic CO bands 
for increasing resolution, from R = 500, which corresponds ap- 
proximately to the value used by McGregor et al. ( 1988) as was 
discussed in Sect.[T| to R = 5000, with the intrinsic CO velocity 
fixed at30kms~'. While for the lowest resolution the 13 COband 
heads are not visible at all, their smooth structure breaks down at 
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Fig. 4. Synthetic, optically thin 12 CO plus I3 CO first-overtone band spectra computed for Tco = 3500 K. The individual panels are 
for increasing spectral resolutions (from top to bottom) and for decreasing 12 C/ 13 C ratios (from left to right). 



the chosen highest resolution, since at high resolution we start to 
resolve individual vibration-rotational lines. The most suitable 
range in resolution that follows from these calculations is thus 
between R = 1000 and R = 3000. 

The terminal velocity of the equatorial winds in B[e] super- 
giants was found to be S> 80 km s (Zickgraf et al. 1 19961 ), so that 
for edge-on systems a resolution as high as R = 5000 would also 
still deliver smooth 13 CO band heads. However, since not much 
is known about the possible inclinations of the galactic B[e] su- 
pergiant candidates, a lower resolution (R 55 3000) is preferable 
for the detection of smooth 13 CO band heads. 



3.3. Influence of optical depth effects 

Besides the CO temperature, the column density of the emitting 
region is a further important parameter that triggers the intensity 
(and thus the observability) of the CO bands. In this section, we 
therefore briefly discuss the influence of the CO column density, 
and hence of the optical depth, on the global appearance of the 
CO bands and, especially, on the appearance of the 13 CO band 
heads. 

We calculate the CO bands for a fixed temperature of Tqo = 
3500 K and carbon isotopic ratio of I2 C/ I3 C = 5 and for increas- 
ing 12 CO column densities. The results for three different values 
of the column density are shown in Fig. [5] Contrary to the mid 
and top panel, which show clear indications for partially opti- 
cally thick emission, the chosen column density in the bottom 



panel is low enough to guarantee that the CO bands are still op- 
tically thin throughout the spectrum. 

The influence of the increasing optical depth becomes obvi- 
ous in the strong increase of the CO "quasi-continuum", start- 
ing especially at long wavelengths, which leads to a reduction 
of the ratio between the peak and the inter-peak (i.e., the quasi- 
continuum) intensity. An additional optical depth effect is the 
broadening of the band heads on their long-wavelength "wing". 
Nevertheless, the 13 CO band heads remain detectable, even for 
significant optical depth values within the 12 CO bands. 



4. Discussion 

The rising of the I3 CO band heads with decreasing carbon iso- 
topic ratio is nicely displayed in Fig. [4] Since the decrease in 
isotopic ratio is a clear indicator of stellar evolution (see Fig. [2] 
and Sect.|3Ji, we can claim that, if 13 CO band head emission is 
identified within the K band spectra of a massive star, it hints 
at an evolved nature of the object. Whether and how this can be 
used to solve the problem of the evolutionary phase of the B[e] 
supergiants, or B[e] stars in general, is discussed in the follow- 
ing. 

4.1. Comparison with non-rotating stellar evolution models 

For a sense of whether the B[e] supergiant candidates in 
Table lA. ll will show 13 CO bands in their K band spectra at all, 
we need to specify those evolutionary phases with suitably low 
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Fig. 5. CO Band structure for increasing CO column density 
(from bottom to top), i.e., increasing optical thickness. 



isotopic ratio and compare them with the locations of the B[e] 
stars within the HR diagram. Such a comparison is performed 
in the left panel of Fig. [6] where we plot the stellar evolutionary 
tracks from Schaller et al. ([1992 ) for non-rotating stars at solar 
metallicity. Those evolutionary phases during which the 12 C/ 13 C 
ratio drops below ^ 20 and further below ^ 15 are indicated. 
Also included are the positions (with error bars) of the galactic 
B[e] supergiant candidate sample. 

Obviously, only the evolutionary tracks for stars with high- 
est initial mass (i.e., M m <: 40 M ) show an overlap of their 
evolutionary phases of low carbon isotopic ratio with the posi- 
tions of the most luminous B[e] supergiant candidates. For these 
targets, the detection of 13 CO bands would immediately clas- 
sify these as B[e] supergiants, and for stars in the mass range 
40 M S= Mji, 5= 60 M Q , even constrain the B[e] supergiant phase 
to late evolutionary phases, i.e., within the second blue super- 
giant phase during the evolution. 

Stellar models with M- m < 15 M reach a carbon isotopic 
ratio of about 15-20 during their evolution along the blue loop. 
Even though there exist some B[e] stars within this mass range, 
the 12 C/ 13 C ratio for these stars might not be low enough for a 
definite 13 CO band detection in low resolution. For these stars, 
a high-resolution high signal-to-noise spectrum covering the re- 
gion around the 2 — > 13 CO band head might be more suitable 
to search for (and resolve) individual vibration-rotation lines. 

For the stellar models of intermediate mass stars (15 M Q ^ 
M in 5 40 M ), the 12 C/ 13 C ratio changes only at the end of their 
evolution, i.e. within the red part of the HR diagram. For this 
mass range thus there exists no overlap of the low carbon iso- 
topic ratio regime with the locations of the B[e] supergiant can- 
didates. 



4.2. The influence of stellar rotation 

It is generally assumed that the disk formation around B[e] su- 
pergiants is linked to rapid rotation of the central stars. In fact, 
for 4 Magellanic Cloud B[e] supergiants, projected stellar rota- 
tion values (i.e., usin/) could be determined. These stars indeed 
seem to rotate at a significant fraction of their critical velocity 
(Gum mersbach et al. [19951 Zickgraf [20001 120061 Kraus et al. 
2008). Whether all B[e] supergiants are rapid rotators, is, how- 
ever, not known and difficult to determine, since for the rest of 
the targets no photospheric absorption lines can be detected due 
to their high-density winds. 

To date, nothing is known about possible rotation velocities 
of the galactic B[e] supergiant candidates. Nevertheless, we in- 
vestigated the influence of stellar rotation on the variation of 
the surface 12 C/ 13 C isotope ratio. For this, we used the tables 
of Meynet & Maeder (2003 ) for rotating stars at solar metal- 
licity. The evolutionary tracks are shown in the right panel of 
Fig. [6] and the evolutionary phases with 20 > 12 C/ 13 C > 15 and 
12 C/ 13 C < 15 are again indicated. The general trend of the stel- 
lar rotation is to enhance the stellar surface abundance with 13 C 
much earlier during the stellar evolution. And interestingly, for 
the rotating models the 12 C/ 13 C isotope ratio drops below a value 
of 20 already during the main-sequence evolution of even the 
stars with lowest (i.e., here 9M ) initial mass, and even stars 
with an initial mass of 25 M experience a decrease in isotope 
ratio to 12 C/ 13 C < 15 during their late evolutionary phase. In the 
presence of stellar rotation, the positions of all our B[e] super- 
giant candidates thus fall into regions with decreased 12 C/ I3 C 
isotope ratio. The search for 13 CO band emission, therefore, 
seems to be a promising way to distinguish between pre- and 
post-main sequence evolution of B[e] stars. 

4.3. The B[e] supergiant candidate GG Car 

In Fig. [6] we emphasized the position of the B[e] supergiant can- 
didate GG Car based on the stellar parameter determinations of 
earlier investigations by McGregor et al. d 19881 ) and Lopes et 
al. ( 119921 ). According to these studies, it seems that GGCar has 
evolved from a progenitor star of M m =* 15 ... 22 M . 

Recently, Domiciano de Souza et al. (in preparation) per- 
formed VLTI/ AMBER observations of GGCar and extracted 
its K band spectrum. Interestingly, the AMBER spectrum of 
GGCar, which has a spectral resolution of R ~ 1500, seems 
to display clearly detectable 13 CO band emission (Chesneau & 
Domiciano de Souza, private com.), classifying GG Car as a B[e] 
supergiant star. Based on the peak/inter-peak intensity ratio, the 
CO emission seems to be rather optically thin, and the detected 
strength of the 13 CO band heads speaks in favour of a 12 C/ 13 C 
isotope ratio of the order of 5= 10. Such a low ratio is, how- 
ever, not expected within the suspected evolutionary phase of 
GG Car, neither from the non-rotating, nor from the rotating stel- 
lar evolution models. This discrepency might indicate that ei- 
ther the luminosity of GG Car (and hence its initial mass) has 
been considerably underestimated, or its progenitor star must 
perhaps have been rotating much faster than 300kms~', result- 
ing in a much stronger I3 C surface abundance enhancement than 
the predictions from the rotating models of Meynet & Maeder 
(2003 ). In any case, the detection of the unexpectedly strong 
I3 CO band heads in GGCar needs to be studied in more detail, 
including a careful modeling of the I2 CO and I3 CO band emis- 
sion (Domiciano de Souza et al., in preparation). In addition, 
more interferometric studies to spatially resolve the CO bands, 
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Fig. 6. Comparison of the positions of the galactic B[e] supergiant candidates from Table lA. ll with evolutionary tracks of non- 
rotating (left panel, from Schaller et al. 119921) and rotating (right panel, from Meynet & Maeder 2003) stellar evolutionary models 
at solar metallicity. The dashed and dotted parts of the tracks indicate those evolutionary phases with 20 > 12 C/ 13 C > 15 and 
12 C/ 13 C < 15, respectively. The position of GGCar is emphasized. 



as has recently been shown by Tatulli et al. (2008), are definitely 
needed. 



5. Conclusions 

For many galactic B[e] stars it is not clear whether they are 
pre-main sequence stars (i.e., Herbig stars) or whether they be- 
long to the group of B[e] supergiants. The main reason for this 
long-standing problem in classification is their often only poorly 
known distance, resulting in huge errorbars for the stellar lumi- 
nosities. In addition, these stars usually have high density cir- 
cumstellar material obscuring the stellar photosphere, so that an 
abundance study of their chemical surface composition is not 
possible. 

To overcome these problems, we suggest a different way 
to distinguish between the pre- and post-main sequence evo- 
lutionary phase of the individual targets based on the increas- 
ing 13 C abundance on the stellar surface at post-main sequence 
evolutionary phases, and hence an increased 13 C abundance in 
the winds of the supergiant stars. Since typical B[e] supergiants 
(such as those in the Magellanic Clouds) show the presence of 
CO band emission arising in their high-density, neutral circum- 
stellar disks, we investigate the conditions under which observ- 
able 13 CO band emission in the K band spectra of B[e] super- 
giant candidates is generated. Only stars for which the 12 C/ 13 C 
isotope ratio on the stellar surface (and hence in the wind mate- 
rial) has dropped substantially, from the initial, interstellar value 
of about 70 to less than 20, are found to be suitable to produce 



observable 13 CO band emission. Such a drop in isotope ratio 
happens naturally in the post-main sequence evolution of ex- 
tremely massive stars, for which mixing processes and the loss 
of the outer shells via strong mass loss result in a strong increase 
in 13 C surface abundance. In the case of rotating stellar mod- 
els, this surface enrichment in the 13 C isotope is much more effi- 
cient so that even stars with moderate initial mass (i.e., stars with 
Mj n ^ 9 M Q ) already show considerable variations, i.e., a drop in 
the 12 C/ 13 C isotope ratio. The computation of synthetic emission 
spectra of both the 12 CO and 13 CO bands for different resolutions 
revealed that a spectral resolution of R ~ 1500 . . . 3000 is ideal 
for the detection of the 13 CO band heads. 



Interestingly, the B[e] supergiant candidate GGCar has re- 
cently been found to show 13 CO band emission, classifying it 
definitely as a B[e] supergiant. The observed strength of its 
13 CO bands is, however, in disagreement with the model predic- 
tions for its postulated position within the HR diagram, meaning 
that either its luminosity and initial mass have been consider- 
ably underestimated, or its progenitor star must have been ro- 
tating much faster than the 300kms _1 used in the rotating stel- 
lar evolution models. And also for the B[e] supergiant candidate 
CPD-57°2874 some tentative 13 CO detection has been reported 
(Chesneau & Domiciano de Souza, private com.), indicating that 
the search for and detection of 13 CO bands is indeed an ideal tool 
to discriminate between the pre- and post-main sequence evolu- 
tionary phase for the still large number of unclassified B[e] stars. 
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Table A.l. Parameters of galactic B[e] supergiant candidates. For each object we provide the ranges in luminosity and temperature 
as found in the literature. We further indicate for each target whether CO band emission has been detected in their K band spectra. 
A question mark means that no CO band emission was detected for R ~ 450. 



Object 


log L/L 


Jeff |KJ 


Ker. 


K band 
spectra 


C(J 


Ket. 


CD-42° 11721 


4.0 + 0.3 
> 4.5 


13-15 000 
16000 


Borges Fernandes et al. (2007) 
McGregor et al. fll988|t 


yes 


7 


McGregor et al. (1988J 


rpn-5?°9?43 


5.9 
5.45 + 0.25 
5.4 


16 200 
13-22 000 
16000 


T nnps et al 1 1 QQ? 1 

J_/U|JCS C L til. I|177Zj|I 

Winkler & Wolf (| 1989) 
McGregor et al. <TT988T> 


yes 


yes 


MrGrepnr et al (19881 

lV±\^VJltll*Jl tl £11. \ 170O| 


CPD-57°2874 


> 4.8 


16000 


McGregor et al. (I988|l 


yes 


yes 


McGregor et al. jl988|l; 




5.7 + 0.4 


17-23 000 


Domiciano de Souza 
et al. \2001\ 






Domiciano de Souza 
et al. (in prep.) 


GGCar 


4.7 


16000 


McGregor et al. (11988b 


yes 


yes 


McGregor et al. (1988); Morris 




5.18 


20 800 


Lopes etal. (1992) 






et al. 11996b: Domiciano de 
Souza et al. (in prep.) 


Hen 3-298 


5.1+0.3 


13-17 500 


Miroshnichenko et al. (2005 ) 


yes 


yes 


Miroshnichenko et al. 12005b 


Hen 3-303 


4.3 + 0.3 


10-20000 


Miroshnichenko et al. (12005b 


yes 


no 


Miroshnichenko et al. (2005 ) 


HD 87643 


4.9 
5.46 
4.0-5.6 


16000 
16 200 
17-18 000 


McGresor et al. 11988b 
Lopes etal. 11992b 
Oudmaijer et al. (|1998b 


yes 


? 


McGregor et al. 11988b 


MWC 84 


< 4.0 
5.4 


18-22 000 
18-30000 


Miroshnichenko et al. (12002b 
Hynes et al. (|2002J 








MWC 137 


4.46 
5.37 


31000 
30000 


Hillenbrand et al. (11992b 
Esteban & Fernandez ( 1998 ) 








MWC 300 


5.7-5.85 
5.1+0.1 


20-22000 
-20000 


Wolf&Stahl 11985b 
Miroshnichenko et al. (12004b 








MWC 314 


5.95 ± 0.2 


16400-29 700 


Miroshnichenko et al. (11998b 








MWC 349 A 


5.7 + 1.0 


20-28 000 


Hofmann et al. 12002b 


yes 


yes 


Geballe & Persson Jl987b; 




6.1+0.1 


32-35 000 


Hartmann et al. 11980b 






Kraus etal. (2000) 


MWC 873 


6.26 
5.0 ± 0.4 
6.0 


13 000 
20000 
11500 


Lopes et al. 11992b 
Miroshnichenko et al. (2003 ) 
Machado et al. (2003b 


yes 


yes 


McGregor et al. 11988b 


W9 (Ara C) 


4.86-5.26 


20-35 000 


Clark etal. 11998b 









Domiciano de Souza, A., Driebe, T., Chesneau, O., et al. 2007, A&A, 464, 81 
Evans, C. J., Crowther, P. A., Fullerton, A. W., & Hillier. D. J. 2004, ApJ, 610, 
1021 

Esteban, C, & Fernandez, M. 1998, MNRAS, 298, 185 

Farrenq, R., Guelachvili, G., Sauval, A. J., Grevesse, N., & Farmer, C. B. 1991, 

J. Mol. Spectrosc, 149, 375 
Geballe, T. R., & Persson, S. E. 1987, ApJ, 312, 297 
Gummersbach, C. A., Zickgraf, F.-J., & Wolf, B. 1995, A&A 302, 409 
Hanson, M. M., Howarth, I. D., & Conti, P. S. 1997, ApJ, 489, 698 
Hartmann, L., Jaffe, D., & Huchra, J. P. 1980, ApJ, 239, 905 
Hillenbrand, L. A., Strom, S. E., Vrba, F. J., & Keene, J. 1992, ApJ, 397, 613 
Hillier, D. J., Lanz, T., Heap, S. R., et al. 2003, ApJ, 588, 1039 
Hofmann, K.-H., Balega, Y., Ikhsanov, N. R., Miroshnichenko, A. S., & Weigelt, 

G. 2002, A&A, 395, 891 
Hynes, R. I., Clark, J. S., Barsukova, E. A., et al. 2002, A&A, 392, 991 
Ignace, R., & Gayley, K. G. 2005, The nature and evolution of disks around hot 

stars (San Francisco: ASP), ASP Conf. Ser., 337 
Kastner, J. H., Buchanan, C. L., Sargent, B., & Forrest, W. J. 2006, ApJ, 638, 

L29 

Kraus, M. 1997, CO band emission from a rotating disk, Diploma thesis, 
University of Bonn 

Kraus, M. 2000, Modeling of the Near IR Emission from the Peculiar B[e] star 

MWC 349, PhD thesis, University of Bonn 
Kraus, M., & Borges Fernandes, M. 2005, in ASP Conf. Ser. 337, The nature 

and evolution of disks around hot stars, ed. R. Ignace & K. G. Gayley (San 

Francisco: ASP), 254 
Kraus, M., & Miroshnichenko, A. S. 2006, Stars with the B[e] Phenomenon, 

(San Francisco: ASP), ASP Conf. Ser. 355 
Kraus, M., Borges Fernandes, M., & de Araiijo, F. X. 2007, A&A, 463, 627 
Kraus, M., Borges Fernandes, M., Andrade Pilling, D., & de Araiijo, F. X. 2006, 

in ASP Conf. Ser. 355, Stars with the B[e] Phenomenon, ed. M. Kraus & A. 

S. Miroshnichenko (San Francisco: ASP), 163 



Appendix A: Galactic B[e] supergiant candidates 

Table lA. ll contains a list of the currently known galactic B[e] 
supergiant candidates. 
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